home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / dscal.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.3 KB  |  61 lines
  1.       SUBROUTINE DSCAL( N, DA, DX, INCX )
  2. *
  3. *     scales a vector by a constant.
  4. *     uses unrolled loops for increment equal to one.
  5. *     jack dongarra, linpack, 3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER           INCX, N
  9.       DOUBLE PRECISION  DA
  10. *     ..
  11. *     .. Array Arguments ..
  12.       DOUBLE PRECISION  DX( 1 )
  13. *     ..
  14. *     .. Local Scalars ..
  15.       INTEGER           I, IX, M, MP1, NINCX
  16. *     ..
  17. *     .. Intrinsic Functions ..
  18.       INTRINSIC         MOD
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       IF( N.LE.0 )
  23.      $   RETURN
  24.       IF( INCX.EQ.1 )
  25.      $   GO TO 20
  26. *
  27. *        code for increment not equal to 1
  28. *
  29.       IX = 1
  30.       IF( INCX.LT.0 )
  31.      $   IX = 1 - ( N-1 )*INCX
  32.       NINCX = IX + ( N-1 )*INCX
  33.       DO 10 I = IX, NINCX, INCX
  34.          DX( I ) = DA*DX( I )
  35.    10 CONTINUE
  36.       RETURN
  37. *
  38. *        code for increment equal to 1
  39. *
  40. *
  41. *        clean-up loop
  42. *
  43.    20 M = MOD( N, 5 )
  44.       IF( M.EQ.0 )
  45.      $   GO TO 40
  46.       DO 30 I = 1, M
  47.          DX( I ) = DA*DX( I )
  48.    30 CONTINUE
  49.       IF( N.LT.5 )
  50.      $   RETURN
  51.    40 MP1 = M + 1
  52.       DO 50 I = MP1, N, 5
  53.          DX( I ) = DA*DX( I )
  54.          DX( I+1 ) = DA*DX( I+1 )
  55.          DX( I+2 ) = DA*DX( I+2 )
  56.          DX( I+3 ) = DA*DX( I+3 )
  57.          DX( I+4 ) = DA*DX( I+4 )
  58.    50 CONTINUE
  59.       RETURN
  60.       END
  61.